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1. Introduction 


The theory of characteristics will he presented generally for 
quasilinear differential equations of the second order in two variables. 
This is necessary because of the manifold requirements to be demanded 
from the theory of characteristics. 

The function x of the two independent variables x, y is assumed 
as satisfying the quasilinear differential equation 


* BX sy * ox yy + 11 - 0 (1) 

This differential equation Is called quasilinear since the highest 
derivatives, namely those of the second order, occur only linearly. The 
coefficients A, B, C, D are functions of x, y, X, x^, Xy The 
differential equations of gas flows set up in chapter I all belong to 


"Character is tihentheorie . " Technische Hochschule Dresden, 
Archiv Nr. 44/2, Chapter II. 
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this type. It trill occasionally he useful to interpret the desired 
integral X(x, y) of the differential equation (l) geometrically as 
Integral surface in the x, y, x - space. 

The characteristics are to he introduced in three ways: 

First, as loci of the possible appearance of small disturbances. 

For gas flows this interpretation is obvious: the characteristic base 
curves here are nothing but the Mach waves obtained in the known 
elementary manner by superposition of sound waves accor ding to Huygens 1 
principle as fronts of a weak disturbance wave. Following it will be 
defined of what type the disturbances or discontinuities are which are 
propagated for instance from the boundary of the region into the 
interior along the characteristics . 

The second interpretation of the characteristics will start from 
the fundamental fact that the characteristics are the sole curves 
from which in general the integral surface x(x, y) can be constructed. 

The third introduction of the characteristics is the one used most 
frequently in mathematical representations. It shows in what sense the 
continuation of an integral surface beyond a characteristic may become 
indefinite. This definition of characteristics is unnecessary for our 
purposes . It is mentioned merely in order to ensure connection with 
the customary mathematical literature; however, omission of this section 
is not detrimental to the understanding of the rest. 

The purpose of this chapter is attained with the development of 
two general approximation methods for the solution of the characteristic 
differential equation system. 


2 . Preliminary Statements 

The characteristics lie on the integral surface X(x, y) . Their 
projections on the x-y plane, or, in other words, their base 
projections, are designated as "characteristic base curves". Let these 
characteristic base curves have the equation 


T](x,y) = Constant 


( 2 ) 


For reasons of a simpler manner of expression, the characteristic base 
curve 


T}(x,y) = 0 


(2a) 
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will "be considered (fig. 1.) The characteristic family of curves (2) may 
he intersected by another family of curves. 


£(x,y) = Constant (3) 


No further data are given concerning this second family of curves; for 
hyperbolic differential equations where two families of characteristics 
appear, the second family of characteristic base curves will be 
selected as I - family. This is mentioned only inc identa lly : for the 
immediately following considerations only t) = Constant are ass ume d as 
characteristic base curves which are intersected by the 
curves I = Constant. Thus one has as coordinate along a characteristic 
curve, I, as transverse coordinate, r\. The derivative of a 
function f(x, y) with respect to g along a characteristic base 
curve (as which q = 0 will be selected below) will be called "interior 
derivative"; to attain it, nothing but the course of the function wi thin 
the considered | - region on T| = 0 is needed. Derivatives with 
respect to t| require knowledge of the behavior of the function to be 
derived outside of the characteristic curve; they will be called 
"exterior derivatives". It is obvious that the conceptions of interior 
and exterior derivatives are very closely connected with the conceptions 
of tangential and normal derivatives. 

After these preliminary statements the announced definitions of 
characteristics are set up. 


3. Characteristics as Loci of Discontinuities 
of the Second Order 

The integral function x and its first derivatives are to remain 
continuous when the characteristic tj = 0 is transversed. Disconti- 
nuities in the second derivatives - the highest ones occurring in the 
differential equation (l) - are to be permissible but with the 
restriction that at least the interior or tangential derivatives still 
remain continuous. The permitted discontinuities concern at most the 
exterior derivatives of the derivatives of the first order. 

These properties are best formulated in the g, q - coordinate 
system. Obviously this is permissible since the geometrical 
interpretations of the conceptions used for definition of characteristics 
are independent of any coordinate system, x is the height of the 
integral surface above the base plane, the first derivatives of x 
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with, respect to any two coordinates define the position- 1 - of the surface 
elements. Interior and exterior derivatives of the first derivatives 
of X, therefore, measure the variation of the position of the surface 
elements along and across the curve t) = 0. 

According to the def ini tion given above X, Xg, X^ are, therefore, 

to remain continuous when rj = 0 is crossed, furthermore the interior 
derivatives of X| and X^, thus Xgg and X^g. The single permitted 

discontinuity of the second order may, therefore, appear only in X^. 

The discontinuity the quantity X^ undergoes in passing from negative 

to positive wild he designated by 

In order to set up the equation of q = 0 in the x, y coordinates 
one has to go back to the original coordinates. 

X xx = + ^Tig^x^x + X gg^x 2 + x r| T lxx + X g^xx, 


X xy - x r)q T lx 1 ly + x 'H|( T lx^y + Hy^x;) + X g|§x£y + Xq^xy + X £&xy> ^ 

x jj = X 4T 1 T iy 2 + 2X qS Vy + x^y 2 + x n T Vy + *s*yy 

Consequently the discontinuities in X^, X^y, Xyy will be 

[xxxj = [ x Tn] nx 2 > 

[XxyJ = [ x qq| ^ 

[ x yy] = [ x r)TiJ % 2 • 


^The position of the surface elements is defined by the direction 
of their normals, the direction- oosines of which in the x-y-X system 
are the proportion X x : Xy : -1. 
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If one now sets up the differential equation (l) for absolutely emal 1 
positive and negative q for the same £ and forms the difference, 
one obtains 

A pxx] + B [ x xyJ + c [ x yy] = °> ( 6 ) 

or according to (5) 

An x 2 + BTj^y + Cq y 2 = 0. (7) 

This equation will be designated as "characteristic condition". 

If one expresses the characteristic base curve q = 0 in the 
parameter form x = x(£); y = y(£), one may write, since the slope 
T l„ 

is — = - for the characteristic condition (7) also: 

X £ V 


Ay 2 - Byx + Cx^ = 0, 


(8) 


where the differentiations with respect to the parameter are denoted by 
dots above them. 

From this quadratic relation one may for hyperbolic differential 
equations where 


B 2 - 4AC > 0 


( 9 ) 


set up the differential equations for two characteristic families of 
base curves. 

Whereas the occurrence of discontinuities of the second order is 
restricted by the characteristic condition, the differential equation (l) 
does not offer a condition for the existence of discontinuities of the 
first order. One can see this readily if one assumes on each side of 
a space curve, selected as carrier of the initial conditions, two 
different positions of the surface elements prescribed for the integral 
surfaces. Discontinuities of -the first order appear in gas dynamics 
as compression shocks. 
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4 . Characteristic Strips as Elements % 

of the Integral Surface 

For the second, interpretation of the characteristics one at first 
orders with regard to a curve tj = 0 on a surface X the positions 
of the surface elements. One then speaks of a strip of the first 
order cc^. The surface x is not a priori assumed to he an integral 
surface. The position of the surface elements is suitably determined 
by the derivatives of X with respect t.o the coordinates | and 
thus X| = p-j_, = q.i* T^ 9 strip is determined in parametric 

representation, thus by x(t), y(t), *(£), <2.^(1 ). 

The problem arises whether it is possible to express the 
quasilinear differential expression on the left side of “the differential 
equation ( 1) : 


AX + BX„, + CX.„ + D 


xx 


*y 


yy 


do) 


merely by the five strip quantities x(g), y(|), x(|), PqU), 1-.CS ) 

and their derivatives with respect to S . One then says that the 
differential expression (10) lies in the strip ct^ or that it is an 

"interior differential expression" pf the strip c^. 

Obviously one may then hope to satisfy the differential equation (l) 
along such a strip without leaving it in a transverse direction, and 
to build up the desired integral surface from such strips. 

According to presupposition, the coefficients A, B, C, D depend 
only on x, y, X, p , q^ which are given directly on the strip. Thus 

it r' emaln s only to be determined when X , X , X can be expressed 

merely by the five strip quantities and their interior derivatives. 
According to the formulas ( 4 ), introducing as far as possible the strip 
quantities, one obtains: 




*xx “ IitVHx 2 + 1 lxl x + £14 ^x 2 + ll^xx + p l^xx 


X xy " + a i4^^y + + p i|^y + ^zy + p l^xy 

V = ‘iinV 2 + 29.16 v 6 y + p iiV 2 + aiT W + Pl| yy 


(11) 
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One sees Imme diately that the sole term which cannot be formed by 
Interior differentiation in <x^ (that is with respect to 5 ) always 

stands first on the right side of these expressions (11). Thus one 
obtains as sole differential contribution of (10) 

1 i n (A \ 2 + B w + °v 2) 

In order to have the differential expression (10) lie completely 
in ct^, it is necessary and sufficient that 

A \ 2 + + Ct ^ 2 = 0 

This is again the characteristic relation ( 7 ). Thus none but the 
characteristic strips can be concerned in the building up of an integral 
surface from strips. Of course. It remains still to be shown that the 
building up from the characteristic strips is actually possible. 


5. Indefinite Continuation of an Integral 
Surface Beyond a Characteristic Strip 

In order to obtain the third customary definition of characteristics 
one starts from a strip of the first order on the integral surface x. 
Let the strip be determined as function of the parameter £ by the five 
quantities x, y, X, X x = p, X— = q.. Due to the meaning of p and q. 
there exists between the five strip quantities x, j, X, p, and q. the 
relation: 


dX _ ciX dx _ 5x dy dr _dy. 

d£ <E d£ + ciy d£ “ ^ 4 d| 

or, if one again denotes the differentiations with respect to the 
parameter by dots above them: 


X = pi + qy 


( 12 ) 


This is the so-called "strip relation of the first order". 

Let the derivatives of the second order which are not given be 
denoted by r = X^, s = X^y, t = *yy The question arises whether 
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it is always possible to determine uniquely with the aid of - the 
differential equation (l) on the strip of the first order the 
derivatives of the second order r, a, and t and furthermore the 
derivatives of higher order so that a continuation of the integral x 
hy a Taylor series appears possible. A well-known theorem of Sonja 
Kowalewskaja deals with this analytical continuation of an integral. 

For the determination of r, s, t from the five strip 
quantities x, y, x, p, q the differential equation (l), for one, 
is at disposal, now written as follows: 


Ar + Bs + Ct = -D 


( 13 ) 


From the meaning of r, s, t one obtains the strip relations of the 
second order: 

d| ~ 5x d| + ciy d£ “ aF aT * 
for which one writes briefly: 

xr + ys = P _ (l4) 


Correspondingly, one obtains as further strip relation: 


xs + yt = q 


( 15 ) 


From the three linear equations 


Ar + Bs + Ct = -D 
ys = p 


xr + 


xs + yt = q 


( 13 ) 

(HO 

( 15 ) 
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r, s, t may be determined uniquely, if the following determinant 
is not. zero: 

ABC 
± j 0 
0 x y 


= Ay 2 - Bsj + Cx 2 


It is interest in g to note that furthermore, in order to make the higher 
derivatives uni quely determinable, non -vanishing of the same determinant 
is necessary and sufficient 2 . 

If the determinant vanishes which thus coincides with the 
characteristic condition derived before (8), r, s, t are no longer 
uniquely determinable, but - if at all - only with exclusion of additive 
solutions of the homogeneous equation system pertaining to ( 13 ), ( l4) , 

(15) . This new criterion of the characteristic condition is, of course, 
very closely connected with the property of the characteristics (used 
in section 3 ) of being geometrical locus for the discontinuities of the 
second order of the integral surface. 

If solutions are to exist at all when the determinant of the equation 
system (13), (l4) , (15) vanishes, additional conditions must exist 
between the coefficients of the left and right sides. Since the rank of 
the matrix 

(ABC - d \ 
x y 0 p 
\0 ± y 4/ 


2 In order to find for instance r z , s z , t^, one may, by 
differentiation of ( 13 ) with respect to x, set up the equation 


Ar x + Bs^. + Ct x 



4* B X S 4" C-g-t 4* 



as well as the two strip relations of the third order 

i*x + ys x = r, 

±s x 4- ytx = s. 


For r, s, t assumed as already determined the condition mentioned 
results for the unique determination of r , s^., t^.. 
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must be smaller than 3, it is necessary that the following determinant 
he zero: 


A C -D 

x 0 p 

0 y q 


= Apy + Cqx + Diy = 0 


(16) 


According to a determinant theorem found hy Kronecker3 this condition is 
sufficient. This will not he discussed in more detail since another 
derivation will he given directly for this second equation (l6) of the 
characteristic strip which is added to the characteristic condition (8). 
As third equation for the characteristic strip one has the strip relation 
of the first order (12) already mentioned: 


X = px + qy 


6. The Characteristic Differential 
Equation System 


The characteristic condition 


Ay^ - Bxy + Cx^ = 0 


( 8 ) 


introduced in three different ways yields the differential equation for 
the characteristic base curves. In order to prepare, according to the 
deliberations in section it-, for the building up of the integral surface 
from strips of the first order on the integral surface x(x, y), one 
has, furthermore, to set up differential equations f.or the height and 
the position of the integral -surface elements which one now visualizes 
as described by X = x(|), X x = p(|), Xy = q(|). Due to the meaning 
of p and q as derivatives of x with respect to x and y, the 
relation 


3 X _ 3X dx 8x d^ 
d| cix d| cjy d£ 


^See for instance M. B8cher: Einfilhrung in die hfihere Alegbra, 

(introduction into the higher algebra), chapter V, 19, first theorem. 
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must be valid, thus, again denoting differentiations with respect to the 
parameter by dots above them. 


X = p± + qj 


( 17 ) 


This is the so-called "strip relation of the first order". 

The given differential equation (l) results after multiplication 
by iy and the designations r = X^, s = X , t = X : 

Arxy + Bsxy + Ctxy + Dxy - 0 (l 8 ) 


According to ; the meaning of p, 4 , r, s, t one obtains the strip 
relation of the second order: 


P = rx + sy, 
4 = sx + ty . 


( 19 ) 


If one inserts the expression for rx and ty from these equations 
into ( 18 ), one obtains: 

(Bsxy - Asy^ - Csx^) + Apy + Cqx + Dxy = 0 (20) 


The bracket vanishes according to the characteristic condition ( 8 ) . 
Thus one obtains as second characteristic differential equation 


Apy + Cqx + Dxy = 0 (21) 


Finally, the third characteristic differential eqmtion is given by the 
strip relation of the first order ( 17 ) : 


X = px + qy 


( 17 ) 
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One therefore has a characteristic differential equation system altogether 
for the five strip quantities x(|), yU), x(|), p(g), q(|): 


Ay 2 - Biy + Cx 2 = 0, 


(22a) 

Apy + Cqx + Bxy = 0, " 

— 

(22b) 

X = px + qy. 


(22c) 


According to presupposition the differential equation (l) is 
hyperbolic, thus 


B 2 - 4AC > .0, 


(9) 


so that the first characteristic condition (8) has two different real 
roots for y: 


y = \ '±, y - A”x. (23) 

One notes for later that A V = £ The parameter of the characteristics 

A 

of the first family is called I, as before, the one on the second 
family is called q. For q - Constant one than obtains just the 
characteristics of the first family, for § = Constant those of the 
second family (see fig. 1). If one substitutes for instance y = A ’i 
into the second characteristic differential equation (21), one obtains 


x(D\ 'x + AA'p + Cq) = 0, 


or 


AX 'p + Cq + Dy = 0. (24) 

The third equation (17) or (22c), respectively, of the characteristic 
differential equation system may remain unchanged. 
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Thus on© obtains for the first family of characteristics the three 
differential equations 


= 0, (25a) 

AX'pg + Cqg + Dyg = 0 , ( 25 h) 

Xg - pxg - qyg = 0 . ( 25 c) 

For the second family of characteristics one obtains, correspondingly: 

y n -X"z T) = 0, (26a) 

AX M PtJ + cq^ + = 0, (26b) 

X^ - PXtj - qy n = 0. (26c) 


If one interprets the differential equations ( 25 ) and (26), 
respectively, as ordinary differential equations for x( £ ) , y ( £ ) , 

X(|), p(|), q(S) and x(q), y(q), X(q), p(ij), q.(rj), respectively, 

one has two indeterminate systems so that from this standpoint one 
cannot arrive at an integration theory of the partial differential 
eq.ua tion (l) . 

However, if one regards the characteristics as Gauss parameter 
curves on the integral surface, so that the latter is described 
by X( £,t]) , x(£,q), y(£,q) one may now interpret the six differential 

eq.ua tions (25) , (26) as partial differential equation system for the 
five quantities x(£,ti), y(|,q), x(|,tj), p(|,T}), g.(|,Tj) which 

seems over determined. 

It is also of importance for the practical calculation that 
the last equation (26c) of the differential equation system is 
automatically fulfilled when the first five are satisfied. It Is, 
therefore, not a case of overdetermination. 

In order to prove this important fact, one multiplies ( 25 b) 
by y^, (26b) by y^ and subtracts: 


A (X 'pgyrj - X"Pf]y§> + C(qgy n - q^yg) = 0. 
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One divides this equation by X’X" = — : 

p l ^ ' *n B- + Vi ' = 0 

Hence one obtains, using (25a) and (26a): 

P I X 71 ■ + Vn - = 0, 

or 

+ sn = + i|y n - (2T) 

If one now differentiates (25c) with respect- to tj: 

% ‘ v* " Vr px ^ - lyiri = 0 

and substitutes the relation (.27) found above: 

x iti - - lyn = °, 

one recognizes that 


- px^ - qy^ 

is a function of tj solely which will be called h(rj). If the strip 
relation (26c) is satisfied on a boundary curve cut from the 
curves tj = Constant (see fig. 2) that is, the first family of 
characteristics, h( tj) vanishes; consequently, the strip relation 
(26c) is satisfied also in the entire interior 1 *-. The condition used 
just now, that the strip relations are satisfied at the boundaries of the 
region, must, of course, be duly'taken into account in prescribing 
initial or boundary conditions. 


^In this conclusion it is assumed that characteristics of the first 
f ami ly which are in the interior of the region somewhere meet the boundary. 
Characteristics the entire course of which lies in the interior of a 
region would represent a rare exception which may be disregarded. 
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The existence of the solutions of the differential equation 
system (25), (2 6) -was first established for a pure initial-value 
problem of (l) by K. Friedrichs and H. Levy, for the mixed problems 
with initial and boundary conditions which are of partic ular import anc e 
for us by F. Frank - .], and B. Aleksejeva. These existence theorems are 
not so interesting at the moment since we may refer to the physical 
evidence of the possibility of solution. Later on, it is true, we shall 
have to deal with the restricting presuppositions made by F. Frank! 
and E. Aleksejeva for their existence proof, since they seem to be 
connected with certain occurrences of physical interest. 


7. Construction of an Approximate Solution of the Characteristic 
Differential Equation System According to the Lattice 
Point and Field Method 

One assumes the six characteristic differential equations of the 
form (25), (2 6) written in differentials: 


dy - X'dx = 0, (28a) 

AX 'dp + Cdq + Ddy = 0, (28b) 

dX - pdx - qdy = 0; (28c) 

dy - X"dx = 0, (29a) 

AX"dp + Cdq. + Ddy = 0, (29b) 

dX - pdx - qdy = 0. (29c) 


The three equations (28) refer to the first family of characteristics 
(slope of the base curves X'), the three equations (29) to the second 
family of characteristics (slope of the base curves X") . One of the 
two equations (28c) or (29c) may be omitted according to the expositions 
of the last paragraph. 

This presentation of the characteristic differential equations 
immediately suggests an approximation method which essentially consists 
in the replacement of the differentials by finite differences. 

One visualizes the x, y - region covered by a net of the two 
families of characteristic base curves. In figure 3 the characteristic 
base curves of the first family with the Blope X' are drawn in solid 
lines, those of the second family with the slope X" in dashed ones. 

The points of Intersection form a point lattice. Accor ding to the first 
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method, the "lattice point method", the position of these lattice 
points and the values of x , p, 1 at the lattice points are to he 
determined approximately. 

A rule for continuation is given by assuming that the procedure 
is performed to include the lattice points 1, 2, whereas the 
continuation is to take place toward the lattice point 3* The indices 
are apportioned according to the numbers of the lattice points. By 
replacing, as announced, the differentials in (28a) and (29a) by 
differences, one obtains for the approximate determination of x, y, 
that is the position of the lattice point 3, the two equations 

73 - yi = x l'U3 - xp), (30a) 

73 - 72 = x 2"( x 3 ~ (3 Ob). 

The approximate values of p^, q.3 at the lattice point/3. are, by 

approximation of the differential equations (28b), ( 29b ) , determined 
correspondingly from the following two equations: 

A l^l ' ( P3 “ Pi) + Ci(q 3 - qi) + Dq(y3 - yq) = 0, (31a) 

"( P3 " P2) + c 2^3 " ^2 ) + d 2^3 " 72) = .°* (31b) 

Finally one may determine Xq according to (28c) or (29c) approximately 
from 


x 3 ' X 1 " Pl( x 3 “ x l) " <11(73 " 7q) = (32a) 

or from 

x 3 - X 2 - P 2 ( z 3 " x 2 ^ " ^( y 3 " 72 ) = 0 (32b) 

The second method, the "field method", offers certain advantages as 
to the representation of the results while the expenditure of calculation 
is the same as in the lattice point method. In the field method one 
starts from the concept that the characteristic base curves divide the 
x,y - region into fields. For every one of these fields an approximate 
value for each p and q is to be determined. The distribution of X 
will then be represented by lines of equally large X (contour lines of 
the integral surface) • According to this interpretation, a calculation 
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scheme Is formed "by the characteristic "base curves; in each compartment or 
field t-wo figures, n am e ly the approximate values of p and q, are written. 

The rule for continuation for the field method is developed from 
that for the lattice point method. One assumes the field method to 
have proceeded so far that the desired approximation values in the 
fields I and H (fig. 4) are already known whereas they are just 
about to be determined for the adjoining field IU. Correspondingly, 
the coordinated lattice point method is assumed to have proceeded to 
include the points 1, 2, 3, whereas it is just about to be applied to 
the lattice points 4 and 5* The data about the approximation values 
obtained from the characteristic differential equations (28b), (29b) 
are of foremost importance; they read, in appropriate sequence: 


A 1 X 1 , (P4 " Pi) + C 1 (q 4 - q x ) + D^yk - y x ) = 0, (33a) 

1 ( Pip - P 2 ) + C 2 (q 5 - q 2 ) + D 2 (y^ - y 2 ) = 0; (33b) 

A 2 X 2 , '(vk ~ P 2 ) + C2^4 “ I 2 ) + d 2(74 " 72 ) = (34a) 

A 3 X 3^ p 5 “ p 3^ + C 3^ a 5 " * 3 ) + D 3^ y 5 " y 3^ = °* (34b) 


By addition of the two equations (33a) and (33b) which have been divided 
by two there follows after appropriate re arrangement 


A-l + Ag X,i' + X.g 1 / Pk + P5 Pi + ?2> 

2 2 ^ 2 2 / 

D 1 + p 2 ( 74 + 7 5 7i + J 2 \ 

+ 2 \ 2 2 / 

A 2 - Aj_ ' + Xg 1 / p^ - P5 Pi " ?2\ 

* 2 2 ^ 2 " 2 j 

A 1 + ^2 ^2* “ X l' /P4 “ P5 Pi " P2\ 

+ 2 2 \ 2 ~ 2 J 

Ag - A x K 2 ' - X-i ' f Pi + p 2 PI 4 . + P5\ 

+ 2 2 V 2 ” 2 / 

. D 2 - ^ ^ 74 " 74 75 - 7 2^ 

+ 2 \ 2 2 /• 


Ci + c 2 f + q.5 q.i + q.2j 
+ 2 \ 2 ‘ 2 J\ 


>(35) 


, C 2 - C 7 ( h ~ ^5 
2 \ 2 


<*5 ~ ^2 j 
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If - In tlie sense of the approximation used. - one now neglects terms 
that are of second, or higher order in the differences, the entire right 
side of (35) is el imina ted and one obtains 


^1 ^2 X 1 ' 


D, 


* \ fo + 
2 V 2 


*- 

L. / ^4 + P5 _ Pi + Pg) + c i + c i /14 + I5 _ gq, + 

V 2 2 7 2 ^ 2 2 / 

v_ + T \ r 


*1 * Is 


9-o- 


(36) 


Wow one identifies certain mean values of the desired quantities in the 
lattice points with the corresponding quantities in the fields to the 
boundary of which the respective lattice points pertain, namely5 


Xi + Xn Xr, + Xq 2Q, + Xc 

- *-g = z i> = *11. - ■ ■ g ^ = *hi; 

2412 . 7l , 412 , 7II , &415 . mi; 


X 1 


2 

+ x 2 


2 

Pi +P £ 


= X.!, 

= PI. 


£ — £ 

Xp + X.o Xh + Xc 

2 ■ x n> i - ” 2 = x ni; 

P 4 **>5 

2 " p ni» 




?2 + ?3 

- PH* 2 


^1 + ^2 ^2 + *3 - % + 

2 = *=*-1' 2 “ 2 = ^-III' 


Finally 


A 1 + ^2 X 1 


+ X 2 * ^1 ^2 *^1 ®2 

. # is transformed 


accord 


'•’The point with the coordinates xj, jj- represents in a certain 
sense the center of the field I, namely the bisecting point of the 
diagonal 1 ... 2. Ihe definition of the field center has to be 
selected for the reason that of -the unknown field HI only the one 
diagonal 4 ... 5, with a course corresponding to that of 1 ... 2, is 
already known. 


♦ • • 
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to identities of the following type: 


A + Ao 

— g-— = A ( X I> 7x> *i) 

A ( x l> Ji, x i, Pi, lq) " Atej, ?i, Xj, Pj, 1j) 
+ g 

72* X 2* ^2* I 2 ) ~ A ( Z I* 7jj ^X> Pj, lj) 


for which one writes abbreviatedly : 


Ai + Aq 
2 


Ai 


A 1 " A I *2. ~ A I 
2 2 


(38) 


If one substitutes this expression and the corresponding ones Into (36) 
and again neglects terms which are of second or higher order in the 
differences, one obtains as a formula of the field method: 


A I^l'(PlII - Pi) + C T (q. TTT " lx) + Dx(7 T tt " 7j) = 0- (39) 


The correspo n di n g formula which can be derived from the equations (34a) 
and (34b) reads: 


A ii x n"(2in " Pn) + c ii(iin " in) + D n(ym " 7n) - 0 (4o) 


It must especially be noted that in (39) the progressing from field I 
to field III takes place by crossing a characteristic base curve of the 
second kind whereas in (39) the slope of the characteristic base 

curve of the first kind appears. A corresponding statement may be made 
regarding equation (40) which regulates the progressing from field II 
to field IH by crossing a characteristic base curve of the first kind. 
For the rest, however, the equations are no more complicated than for 
the lattice point method. 

The closing of the new field III by characteristic base curves 
connecting the lattice points 4 and 5 with point 6 would be caused 
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according to the lattice point method by 




yg " y 4 = V( x 6 " x 4) 

76 - y 5 = ^5 ,r (x6 - x 5 ) 

If one replaces in the field method by X TTT ', X^ M by X TTT ": 


(41) 


y 6 ' 


y 6 - 


74 - X m '( x 6 " x 4 ) 


75 - ^TTT ( x 6 x 5 ) > 



(42) 


one has altered the equations (4l) only by quantities -which are of 
second order in the differences, so that the equations (42) can be used 
as further equations of the field method. 


The contour lines of the integral surface, that Is, the lines of 
equally large x, are constructed simply according to the relation 


dx = pdx + qdy = 0 


(43) 


as lines of the slope 


_ E 

dx q 


(44) 


starting from an initial distribution of X. In every field p and q 
remain unchanged according to the field method so that -one obtains an 
approximation of the contour lines of the integral surface by series 
of lines. 

The equations (39)> (40), (42), and (44) determine the field 
method. They are no more cumbersome than the equations of the lattice 
point method. 

So far the progressing had been represented only in the interior 
of the region. How the two methods have to be altered for the 
boundary will be indicated by an example. 
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Let the lattice point method "be performed Including the points 
1 , 2 , 3 and "be just about to he extended to the next lattice point 4 
(point of intersection of a characteristic base curve of the second 
kind with the heavily drawn boundary, figure 5). According to (29a) 
one then obtains approximately 


74 " ^3 = x 3"( z 4 “ *3)- 


The second equation for determination of the coordinates x^, y^ of 

the lattice point 4 is yielded by the equation of the boundary curve. 
According to (29b) one obtains further approximately: 


A 3 X. 3 "(P 4 - P 3 ) + C 3 (q 4 - q 3 ) + D-^y^ - y 3 ) = 0. 


The second equation for the determination of p^, q^ must be given by 

the boundary condition. X 4 may then be approximately determined, 
according to ( 29 c), from 


X 4 - x 3 - p 3 (x^ - x 3 ) - q 3 (yh - y 3 ) = 0 . 


According to the field method one lays a characteristic base curve 
of the second kind through the center of the field I (bisecting point 
of the diagonal 1 ... 2 ) to the point of intersection with the boundary. 
This point of intersection has been marked in figure 5 hy the field 
number II. Within the accuracy of the field method this characteristic 
base curve runs parallel to 2 ... 3 • Then one obtains according to 
(29b): 


A lV' (p H ' P I^ + °l^n " V + ^II ' *j) = °* 

The second relation necessary for determination of Pjj, Ijy must 

be yielded by the boundary condition. The closing of the field II 
by a characteristic base curve of the second kind and the plotting 
of the curves of equal X - values takes place in the customary manner . 

The methods described yield only a first approximation. Tor lin ear 
problems improved approximations are often obtained by a refined approxi- 
mation of the differentials in the characteristic differential equation 
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system (28), (29) with the aid of higher differences. For nonlinear 
problems such an attempt does not seem very promising. Iterative methods, 
on the other hand, will probably lead, with tolerable expenditure of 
calculation, to Improved approximations and finally also to an estimate 
of errors. 


8. Literature 

1. The development of the theory of characteristics given in 
sections 1 to 6 is governed entirely by the needs of the practice and 
thus deviates from the usual text book representations. However, 
similar interpretations may be found represented particularly in the 
book by R. Courant and D. Hilbert entitled "Methoden der Mathematischen 
Physik II" (Methods of Mathematical Physics II), Berlin 1937. Hence, 
for instance the concepts of interior and exterior differentiation 
have been taken over. 

2 . The existence proof for the pure initial problem has been 
given by K. Friedrichs and H. Lewy in Mathematische Annalen, Bd. 99, 
p. 200, 1928. Accounts of this existence proof may be found, for 
instance, in J. Eadamard's "Lecons sur le probldme de Cauchy," (Lessons 
regarding the problem of Cauchy) p. 487, Paris 1932 and in Courant - 
Hilbert l.c., p. 326. For the mixed problems occurring in the 
applications in gas dynamics the existence proof lias been given along 
the same lines by F. Frankl and R. Aleksejeva in a report entitled 
’^Two boundary -value problems from the theory of the hyperbolic partial 
differential equations of the second order with application to gas 
flows at supersonic velocity, " (Russian) , Matematiceski Sboraik, 

issue 1934, p. 483. 

3. The distinction between lattice point and field method has been 
newly introduced here. A report on the lattice point method which was 
developed particularly by I. Massau, Gent 190G. to 1903, is to be 

found in the encyclopedia article by C. Runge and Dr. A. Willers II C 2 on 
numerical and graphic integration, p. 160. The field method had so 
far been developed only for a special case by L. Prandtl and A. Busemann, 
"W&Lerungsverfahren zur zeichnerischen Ermittelung von ebenen Strimungen 
mlt Uberschallgeschwindigkeit, " Stodolaf estschrift Zurich 1929 
(Approximation method for graphic determination of two-dimensional flows 
of supersonic velocity, Stodola anniversary publication Z’Hrich 1929) . 

An essential simplification in this special case is based on the fact 
that the differential-equation coefficients A, B, C depend only on p 
and q whereas D altogether vanishes. For the general case treated 
in this report there had to be found, in contrast, above all a significant 
definition of the field center. 
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Append It 

Simplified Derivation, of the Field Method 


Whereas the field method ■was developed from the lattice point 
method in section 7 of chapter II, it is to "be derived here directly 
from the characteristic differential equations. 

According to the "basic concept of the field method the i , y - 
region Is divided into fields In the manner of a calculation scheme 
"by a net of characteristic "base curves^ the appr oxima te values of p 
and q. Eire written in. these fields. The method is assumed to have 
proceeded to a point where the approximate values in the fields I and II 
(see fig. 6) are already known whereas they are just about to "be 
determined for the adjoining field HI. The desired rule for continuation 
will be obtained by means of an appropriate definition of the field 
centers . 

Since of the unknown field HI 4- ... 5 is known as the sole diagonal, 
the bisecting point of this diagonal Is selected as field center. The 
field centers of the fields I and H are then defined as bisecting 
points of the correspondingly situated diagonals 1 ... 2 anfl 2 ... 3 . 

One assumes the field numbers I, II, -III written at these field centers. 

The comer points 1 ... 4- are, within the scope of our approximation, 
connected by a straight line. Since this 1 ire is clos ing the field I, 

Its slope is selected to equal that of a characteristic base curve 
corresponding to the approximation values prevailing in the field I, thin : 

74- “ 7i = " x i)- (4-5a) 


Similarly: 


7^ 7g — ^*tt ' “ "^ 2 ) • 


(^5b) 


According to the definition of the field centers there is 


Xi + Xo 

*1 = — 5 — > yi 


71 + 72 


2 

x^ + x, 


Xtt — 


y A n 


* iL + x 3 y 
2 ' 


72 + 7o 

— 5“^* 


L IH 


5 J 4 + 

» ym = — 2- 


> (46) 
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Adding (45a) and (45b) and dividing by 2, one obtains 


y 4 + y 5 *1 + J 2 , t x 4 " X 1 , , x 5 - 3Cg 

2 " 2 “ X I 2 ^ X H 2 

2 2 / 


— Xn 


x 4 + 


(^u' - - X J 


U47) 


Neglecting terms which are of the second order in the differences one 
obtains according to (46) an a (47): 


y TTT " y l “ X l'^ x III ' Z I^* 


(48a) 


Correspondingly, it can be shown that 


7 ttt ■ 7n = x n"( z in “ x n) • 


(48b) 


Therewith it is proved that the lines connecting the field centers in 
first approximation are characteristic base curves. 


For the transition from I to III or II to III, respectively. It Is 
now permissible to set up the characteristic differential equation ( 28 b) 
or ( 29 b), respectively, with the differentials replaced by differences. 


One obtains as rule for continuation for the determination of 


and q, 


^III 


III' 


A t X t '(P ttt " Pi) + C T (q TTT - qj) + DjCyjjj - yj) = 0, 

a ii x ii^piii “ ph) + c n^in " ^n) + D n^ y m " ^n) = °- 


> (49) 


Field III Is then approximately bounded by characteristic base curves: 


76 ~ 7k = ' ( x 6 - x 4^ 

y 6 “ y 5 = X Hl " (2: 6 " x 5 } > 



NACA TM 1242 


25 


if 6 is the ind ex of the unknown corner point of field III. The 
corresponding rule had heen assumed "before for fields I and II. 

The construction of the contour lines (X = Constant) of the 
integral surface has already "been discussed on page 20 (equations (43) 
and (44)), the modifications of the method becoming necessary at the 
boundary on page 21. The representation of these facts is not affected 
by the new viewpoint. 


Translation by Mary L. Mahler, 
National Advisory Committee 
for Aeronautics. 
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Figure 4 



Figure 5 



Figure 6 


NACA - Langley Field, Va. 



